Fractional dynamical systems defined on 
fractional jet bundle and applications in 

economics 

■ Mihai Boleantu 

O 



o 



Dept. of Economics, West University of Timisoara 
I email: mihai.boleantu@fse.uvt.ro and 

miliaiboleantu01@yalioo.com 

O 

CN . Abstract 

Using Caputo fractional derivative of order a we build the fractional 
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^ ! 1 Introduction 

, Fractional integration and fractional differentiation are generalizations of the 

' notions of integer-order integration and differentiation, and include n-th deriva- 

, fives and n-fold integrals as a particular case {n denotes an integer number). 

■ Fractional calculus has been used successfully in various fields of science and 

engineering. These applications include classical and quantum mechanics, field 
theory, optimal control and economics. In the past few decades many authors 
[1] ) [H] ) |10| have pointed out that fractional-order models are more appropriate 
than integer-order models for the description of memory and hereditary prop- 
erties of various processes. This is the main advantage of fractional derivatives 
}J] ' in comparison with classical integer-order models in which such effects are, in 

5^ \ fact, neglected. For example, the advantages of fractional derivatives become 

apparent in modeling mechanical and electrical properties of real materials, as 
well as in the description of rheological properties of rocks (see [2^, '^). Also, 
in economics, fractional calculus has been used for deriving an analytical model 
of the tax version of the Fisher equation that incorporates a memory function 
for stock prices and inflation rates (see [1]). 

The problems have been formulated mostly in terms of two types of frac- 
tional derivatives, namely Ricmann-Liouvillc (RL) and Caputo [2], [1]. Among 
mathematicians, RL fractional derivatives have been popular largely because 
they are amenable to many mathematical manipulations. However, the RL 
derivative of a constant is not zero, and in many applications it requires frac- 
tional initial conditions which are generally not specified. In contrast, Caputo 
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derivative of a constant is zero, and a fractional differential equation defined in 
terms of Caputo derivatives require standard boundary conditions. For these 
reasons, Caputo fractional derivatives have been popular among engineers and 
scientists. 

A simple and really geometric interpretation of several types of fractional- 
order integration is given in iQj. Based on this, a physical interpretation of the 
Riemann-Liouville fractional integration is proposed in terms of inhomogeneous 
and changing (nonstatic, dynamic) time scale, and it is shown that the Caputo 
fractional derivative has the same physical interpretation as the Riemann- Li- 
ouville fractional derivative. 

In this paper, using Caputo fractional derivative, the fractional jet fibre bun- 
dle is built on a differentiable manifold and its main geometric structures are 
emphasized. Some significant examples from economics are presented. In sec- 
tion 2, some useful properties of the Caputo fractional derivative are recalled 
and the fractional osculator bundle of order k is described. In section 3, the 
fractional jet fibre bundle is defined and the fractional Euler-Lagrange equations 
together with the fractional Hamilton equations are established. In section 4 
two fractional economic models are studied. 



2 The fractional osculator bundle of order k on 
a differentiable manifold 

2.1 Caputo fractional derivative of order a 

Let X : [a, 6] ^ M be a smooth function, a e M, a > 0, and m G N* such that 
m — 1 < a < m. The left (right) Caputo fractional derivative of a; ^24 is the 
function 

a^t -^It; — r(m-Q) Ja (t-s)° + i-™"* 

.D°'x(t) - 1 f ^ (-D-.'-'C.) ■ ^ > 

where x^™^(s) — j-^x{s) and T is the gamma function of Euler. The following 
properties result from ^ (see [3^). 

1. aD^{ciXiit)+C2X2{t)) = CiaDfXi{t)+C2aDfX2{t), Vci,C2 G M, 

Vxi, X2 : [a, h] —> M. 

2. D],Dfx{t) - aD^+^x{t) + |^^x(™)(a), D] = f^. 

3. If {a„} 

n>o ^ ^^^^ number sequence with lim — p G N*, then 

— n — >-oo 

lim {aD'^-x){t)^Dlx{t). 

n^oo 

4. (a) If x{t) ^c,yte [a, 5], c G R, then aDfx{t) = 0. 
(b) lix{t) = t^, Vi G [a,b], ,Dfx{t) = 

5. If xi^X2 ■ [a, 6] M are analytic functions, then 

,A"(^ia;2)(t) = E u KA""'a;i(i)A'=X2(t), 



fc=0 



where A' = Tt°-°Tf 



2 



6. {aDfX2it))dt = -/>2(i) iaDfx,{t))dt. 

7. If a; : [a, b] —> R is an analytic function and G (a, 6) then 



2.2 Higher-order fractional osculator bundle 

Let a G (0; 1) be fixed and M a differentiable manifold of dimension n. Two 
cm^ves ci,C2 : / — > M with Ci(0) = C2(0) = G Af, G /, have a fractional 
contact a of order fc G N* in xq, if for any / G J-{U), xq E U , U a. chart on M, 
it holds that 



A""(/oci)I*=o =i?r(/°C2)|t=o , (3) 

where a = and Df' = o-Df". 

The set of equivalence classes defined by ([3]) is called the k- osculator frac- 
tional space at M in xq and it will be denoted by Osc"^{M). If the curve 
c : / — > Af is given by x'' — x'^{t), t £ I, i = 1, n, in the chart U, then the class 
[clS^'^ G 0<„'=(M) isgivenby 

a;«(t) = a:* + ^^^^rV^"-), tG (-£,£), 

where i = a = l,fc and a;* = a;*(0). The fractional osculator bundle 
of order k is the bundle {Osc°'''{M), M) where Osc"'^{M) = U^gm Osc^''(M) 
and tt"''^ : Osc°'^{M) ^ M is defined by 7r"'=([c];^'=) = x, V[c]^'= G Osc°'^{M). 

For / G J^{U), the fractional derivative of order a, a G (0, 1), with respect 
to the variable is defined by 



fU) ^ r dfix\...,x^-\s,x^+\...,x-) 1 

where x G C/ab — {x £ U, a* < < 6', « = 1,"}, x^, i = l,n, are the 
coordinate functions on U and {^}, * = 1,", is the canonical base of the 
vector fields on U. 

Using the fractional exterior differential [IJ d" : J-{U) — > V^iU) given by 

d'' ^d{x^YD'^,, (5) 

where (a;-' )" G T{U) and 'D-^{U) is the module of the differential 1-forms on 
U with the canonical base {da;*}, i = l,n, we get [3] 

Proposition 1. (a) With respect to the transformation of coordinates = 
a;'(a;^, a;"), i — l,n, det ^ffj^ 7^ 0, corresponding to the charts U, U' , U n 
U' ^ 0, we /lave </ie relations 

d{xT = J}{x,x)d{xir, D^. = ji{x,x)D^,, (6) 
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where J^^ix,x)^T^D^,{x=r. 

(b) The transformation of coordinates on (7r"*')~^(?7 n f/') C Osc"''{M) are 
given by 

r{a{a - ^ ^ 



r(a) 

r(a) ^ r(a) 



3 The fractional jet bundle on a differentiable 
manifold. Geometrical objects 

3.1 The fractional jet bundle 

The fractional jet bundle of order a on the manifold M is the space J" (R, M) = 
R X Osc°'{M), where Osc°'{M) is the fractional osculator bundle. The triplet 
(J"(R, M), ttq, M) has a structure of differentiable fiber bundle, where ttq : 
J"(R, M) — » M is the canonical projection. If (a;'), i = 1, n, are the coordinate 
functions on the chart U C M, then the coordinate functions on {t^q)~^{U) are 
given by {t, x\ y'^")) where = Y(T+^Dfx'{t) \t=o , i = T~n. 

From the properties of the Caputo fractional derivative (subsection 2.1) and 
from Proposition 1, it follows: 

Proposition 2. Let us consider the functions (t)", (a;*)", of 
Hi^SrHU)), the 1-forms r^d{tr , Tirr^d{xT, rpT^c^((y'("))") of 
X>i((7rg^)-i((7)) and the operators Df, D'^i, D^^t^), i = l,n. The following 
relations hold: 

r^d(P)(z??) = 1, ^^d{xnD'^^,) = 5i, 

^^d{y^M)»{D»,^^)=S^j. 

The module generated by the operators D^, D"i, D'^^^-,, i = l,n, will be 
denoted by X" {{Trg)~^ {U)) . For a ^ 1 this module represents the module of 
the vector fields defined on tTq ^(i7). 

Let us consider two charts U, U' on M with ?7 n C/' ^ 0, {tt§)-'^{U), 
{■jT^y^iU') C J"(IR, M) the corresponding charts on J"(R, M) and the coor- 
dinate functions (x*), (x^), respectively, (t, y*^"^), {t, x^ ,y^^"^). From Propo- 
sition 1, we obtain the transformations of coordinates 

X* — x'^(x^, x") 

(8) 

t/*(") = JUx,x)y^^"\ 
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3.2 Geometrical objects on J"(M, M) 

On the manifold J"(R, M), the following canonical structures may be defined: 

1 

§ = 9' 

From ([H]) and from Proposition 1, it follows that the structures ^ have 
geometrical character. 

The vector field f e X"{{n^y^{U)) is called a /rarfionaZ vector field (FVF) 

iff 

d(i)"(f) = l, 0\T) = O, i = T;7i. (10) 
In local coordinates, (FVF) is given by 

T = D^+y<'^^D^,+FWp,^,, (11) 

where F' e C°°(«)"i(C/)), i = 1771. The integral curves of {FVF) are the 
solutions of the fractional differential equations (FDE) 

D^°'x\t) = F\t, x{t), D^x{t)), i = T~fi. (12) 

The system (jl2p . with given initial conditions, admits solution [S]. 

Let L £ C°°(J"(R, AI)) be a fractional Lagrange function. By definition, 

the Cartan fractional I- form is the 1-form 9l given by 

o'l = Ld{tr + SiL). (13) 
The Cartan 2-form ljl is defined by 

L^L^d^H, (14) 
where d" is the fractional exterior differential 

= ^(t)"^" + d{x')°'D'^. + d(y*("))"£l^.(„,. (15) 

In the chart (ttq )^^(?7), 6'l and lj^ are given by 



dL^{L^ ^^^fi-)D-^,^,{L))d{tr + ^^^D-^^^,{L)d{xT 

a?i = ^^d(t)" A d(a;*)" + S,d(t)" A d(y'("))"+ (16) 

Ayrf(a;*)" A d(a;J)" + B,jd{x^)" A %J("))", 

where 
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Proposition 3. // the fractional Lagrange function is regular i.e., 

det {jhj^ 0, gij = D'^^^^^D'^^^^^^L, then there is a fractional vector field (FVF) 

such that i Q (oji) = 0. In the chart {t:^)^^{U), Tl is given by 

Tl = + + W D^,,^, , (17) 



where 



AP=g^'^{D^,L-dfiD^,,^,L)) 



(18) 



3.3 The fractional Euler-Lagrange equations 

Let c : t e [0, 1] ^ 6 Af be a parameterized curve such that Imc C 

?7 C M. The extension of the curve c to J"(R, Af) is the curve c" : t ^ 
[0,1] (t,x'{t),y'(°'\t)) e J"(R, Af) with /mc" C (vr^)-^^^) C J"(R, Af). 
Let L G C°°(J"(M, Af)) be a fractional Lagrange function. The action of L 
along the curve c" is 

A{c")= [\{t,x{t),y^^Ht))dt- (19) 



Let Ce : i e [0, 1] (x'(t,e)) G Af be a family of curves, with e sufficiently 
small in absolute value so that Imce C Af , co(i) = c{t), f?"ce(0) = f?"ce(l) = 0. 
The action of L on the curves c" is 

^(c^)= r L{t,x{t,e),y("\t,e))dt, (20) 



where (t,e) = r(i+a) -^"^'ft' * ~ '^''^'^ action (1^ has a frac- 

tional extremal value if 



D^A{cf)U^o =0. (21) 

Using the properties of the Caputo fractional derivative (subsection 2.1), it 
results 

Proposition 4. (a) If the action (j20p reaches a fractional extremal value then a 
necessary condition is that c{t) satisfies the fractional Euler-Lagrange equations 



^ + y'(")f?;;^. + y'(^"^f?" („) 
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(b) // the fractional Lagrange function is nondegenerated, then the equations 
(j22p are the fractional differential equations associated to the fractional vector 

a 

' ' Tl given by (|17p . 

(c) // the fractional Lagrange function is nondegenerated, then the system 
P2|) may he written in the form of the fractional Hamilton equations 

D^p'f^ = -D^.H, D^x^ = D^^^,H, (23) 

where 

H=p'^^D^x^-L{t,x{t),y('^\t)) ^^^^ 
P^^ ^Dp^^,Lit,x,y(^^), z=T;7I. 

(d) If for /, h : J"'^(R, M)* — > R the fractional Poisson bracket is defined by 

{/, hr = D^^.^JD^.g - D^,fD^^,^,g, (25) 
where the local coordinates on J^(M., M)* are {x,p^°''^), then 

{H,p^^r^Dy'f\ {H,xT^D^^\ ^^h^- (26) 

4 Economic models described by fractional dif- 
ferential equations 

4.1 The fractional model of Liviatan-Samuelson 

Let us consider the fractional Lagrange function L G JF(J"(M, Af)) given by 

L{t,x,y^^^) = Li(x,2/("))ii;„(-pr), (27) 

where Ea is the Mittag-Lefflcr function, Ea(t) — J2T=q r(i+ak) 
is the discount rate. Using the relation DfEai—pt"") = —pEa{—pt'^) and 
Proposition 4 we obtain 

Proposition 5. (a) The fractional Euler- Lagrange equations (|22p for (|27p are 

y^'-^'^^Dp,^,Dp,^,L, + - pDp,^,L, - D^.L, =0, z = 1771. 

(28) 

(b) IfLi e T{J"{R, R)) is of the form 

Li(x,y(")) = t/(.g(a;)-y(")), (29) 

where U is the utility (welfare) function and c — g{x) ~ j/^"-* is the consump- 
tion function, then the fractional Euler- Lagrange equation is 

r(l + a)^U"{g{x) - y("))(y("))2a _ r(l + a)D'^g{x)U"{g{x) - y("))y(") 
+pU'{gix) - y("))r(l + a) - U'{g{x) - y^^^)D^g{x) = 0. 

(30) 
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Proposition 6. If Li G J^(J"(M, R)) is given by 

ii(x,y(")) = -ai(2/("))2"-a2(2/("')"a;"-a3x2", 01,02,03 gR, (31) 
then the fractional Euler- Lagrange equation is 

air(l + a)r(l + 2a)?;(2") _ (a2r(l + a)^ + pair(l + 2a))(j/("))" ^^^^ 
+a2r(l + - (a3r(l + 2a) + pa2T{l + a)^)^" = 0. 

If ai = 03 = ^, 02 = a and a — > 1, then the function Li and the equation 
([32|) become 

Li{x, i) — ~ ~ ~ 2^^' ^^"^"^ 
The equation represents the classic model of Samuelson [7]. 

4.2 Fractional economic models with restrictions 

Let us consider the Lagrange function L e C°°(J"(M, M)) and the function 
F G C°°(J"(R, M)). The fractional Euler-Lagrange equations of L on the 
restriction F(a;,y'^"^) — 0, G (ttq )^-'^([/), are given by the fractional 

Euler-Lagrange equations of the fractional Lagrange function 

L2{t, A, X, = L{t, X, y^^y) + XF{x, y^")), (34) 
where X(t) is a Lagrange multiplier. From Proposition 4, we obtain 
Proposition 7. (a) T/ie fractional Euler-Lagrange equations of p4p are; 

D^.L + XD-F - dr{D-^,^,L) - Xy^(-)D-,{D-^,^,F) ^^^^ 

(b) // the Lagrange function is given by (j27p t/ien t/ie fractional Euler- 
Lagrange equations (155)1 become 

E^{^pt-){D-M-^pDp,^,L, - y^('^)D-,{Dp,^,L,) 
-2/^('"'i^;^.(„,(^;;.(c,ii)) + A(i?^.F - y^-(")i^^,(i^;^,,„,F) (36) 
-y^(2a)^a^^^pa^^^^)) _ D^XD^^^^.F = 0, 

/or i = l,n. 

The fractional model of investments with restriction is described by the 
function Li{K,I,N) where K{t) = x'^{t), I{t) = x'^{t), N{t) = a;3(i) represent 
the capital stock, the investment and the labor, respectively. The restriction is 
given by F{K'^°'\ K, I, N) = (t>{K,I,N) - = 0, = From jSH) 

we obtain the fractional Euler-Lagrange equations 

E^{~pt'')D'^Li + XD'}^<j> = -DfX 

Ea,{-pt")DfLi + XD%(l)^0 (37) 

Ec,{^pt'')D'^Li + XD%(l>^0. 
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If Li and (j) satisfy the relations 

K'^D$,L,+I-DfL,+N-D-^L, = ^^rL,, r eW ^^^^ 

from ([57]) we get 

rEai~pt°')Li = -AiiT^") - r(l + a)K"X^°'l (39) 
For a ^ 1 the classic model of investments [7 is obtained. 
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